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ABSTRACT 


The Thermal Analysis Variant of the COMOC (Computational Continuum Mechanics) com- 
puter system solves problems involving transient heat conduction and convection in stationary con- 
tinua spanning arbitrarily irregular two-dimensional and axisymmetric solution domains. COMOC is 
based upon a finite element solution algorithm for the energy equation, and solves for the transient 
nodal temperature distribution using a highly stable and automatic explicit integration procedure. 
COMOC is extensively user-oriented, requires minimal input, and no a priori knowledge concerning 
the stability character of the differential equation system. It can readily output computed data in 
user-specified format fields, that geometrically resemble the solution domain discretization (for rapid 
engineering evaluation). This report provides the user with complete information for applying 
COMOC to a specific problem. » - 
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1. SUMMARY AND GUIDELINES 


The Thermal Analysis Variant of the COMOC (Computational Continuum Mechanics) computer 
system can solve diverse problems involving transient heat conduction and transfer in stationary con- 
tinua spanning arbitrarily irregular two-dimensional and axisymmetric spaces. COMOC embodies a 
finite element algorithm for solution of the energy conservation equation, and determines the transient 
nodal temperature distribution using a highly stable explicit integration algorithm. COMOC is highly 
user-oriented, including both input and output phases, and solution accuracy is automatically main- 
tained while minimizing computer execution time. Various manual modes of operation of COMOC 
are user-specifiable through parameter overrides, to adapt the code to any specific requirements of 
a given problem. A complete description of the theoretical basis for COMOC, as well as detailed 
discussion concerning the various technical features of the program, are presented in Bell Report 
9500-920257.* 

The following general guidelines will assist the user in adapting COMOC to a specific 
problem. 

A. DISCRETIZATION 

Up to 1000 finite elements, in two-dimensional or axisymmetric space, involving up to 700 
node locations may be specified. The solution domain closure and the discretization may be ar- 
bitrarily irregular. The nodes of the triangular shape finite elements must sequence in a counter- 
clockwise fashion. 

B. PROBLEM DESCRIPTION INPUT 

Input phases occur at the node and element levels. Any specific initial temperature distri- 
bution is accepted, and up to 100 nodes may have fixed temperatures (identified by inputting the 
negative of the desired temperature). Internal heat generation is specifiable in every finite element. 
Each element may have one side (only), connecting the first two nodes listed for that element, on 
the global domain closure, where any combination of convection or fixed heat flux boundary con- 
ditions may be applied. In two-dimensional problems, convection or heat flux boundary conditions 
may be applied to the element surface as well, and each element may have a unique thickness. 

C. THERMOPHYSICAL PROPERTIES 

The thermal conductivity may have tensor properties, and all thermophysical data are 
specifiably temperature dependent via tables with up to twenty entries. In addition, the solution 
domain may consist of up to six different materials, each with its own thermophysical specifi- 
cations. Each finite element of the discretization must be identified with only one material 
property description. 


*Baker, A. J. and P. D. Manhardt, “Finite Element Solution for Energy Conservation Using a 
Highly Stable Explicit Integration Algorithm,” Bell Aerospace Research Report 9500-920257, 
October, 1972. 



D. INTEGRATION INPUT 

Only the initial and final values of the integration interval, the integration option (1 or 2) 
in COMOC and the desired output print interval need be specified. COMOC will automa tica y 
ex trele i„.e g ra,,on step size wh.le maintaining a stable so, u, ion. Other input - ures a low 
various levels of manual control, including tighter accuracy control, marching at a fixed step , 
specifying an initial minimum step size, and output at exact time stations. 


E. 


OUTPUT 


The standard subroutine provides tabular output of temperature and time derivative of 
temperature at the nodes for the specified output interval. The user may readdy prepare a 
special subroutine to output these variables in a format geometrically similar to the Physical 
shape of the problem. In either case, output is typically provided at intervals close to the 
quested value. By setting a parameter, and at the expense of slightly longer computer run. ti , 
output is provided at exactly the specified intervals. 

Two additional intermediate output options exist. By setting a parameter complete 
printout of all matrices computed for each finite element is obtained. Additional y, comp e e 
printout of the matrix operations and assembly of the global solution vector, associated with 
the derivative evaluation operation within the integration algorithm, may be obtained for a user- 
specified number of passes through COMOC. These outputs are useful for checking that pro- 
Mem input has been properly prepared. All the output features are illustrated in the sample 

problem discussed in this manual. 

COMOC employs non-dimensional dependent and independent variables throughout-, therefore, 
input in any consistent units system is acceptable, and output will occur in that same system e 
Hollerith statements in COMOC assume use of English units (Btu, hr, ft R). A change to IS , 
operation in a dual output mode (English and ISU), could be readily made if required. 
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II. DESCRIPTION OF ROUTINES AND FLOW CHARTS 

A macro flow chart illustrates the overall flow of the computer program. This is followed 
by descriptions and flow charts of all the subroutines which constitute the Thermal Analysis 

Variant of COM OC. 
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MAIN 

This is the executive routine which performs all input and nondimensionalization opera- 
tions and initiates the execution of the geometry and integration routines. 


5 



MICRO FLOW-CHARTS 


READ TITLES 
AND 

GENERAL 

INFORMATION 


CALL 

GEOMETRY 

PACKAGE 


READ 

ELEMENT 

INFORMATION 


NONDIMENSIONALIZE 

GEOMETRY 

MATRICES 


READ 

NODAL 

INFORMATION 


READ ARRAY 
TO 

ORDER 

OUTPUT 


TEMP < 0.0 ? 


TEMP = 

-TEMP 




STORE 

NODE NO. AND FIXED 
TEMPERATURE 
VALUE 


READ 
MATERIAL 
CP, RHO 
AND K TABLES 


CALL QKNINT 
TO 

BEGIN TIME 
LOOP 




COMPUTE 

NONDIMENSION ALIZING 
FACTORS 


LAST 
CASE ? 


READ 

INTEGRATION 

INPUT 


6 











DERIV 


This routine computes and assembles the derivative of temperature which is required by the 
integration package to evaluate the temperature. 

It is called from QKNINT and QKNUIN. 
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FECOM 


This routine prints the title page for each case. 




It is called from MAIN. 
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FECOM 
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GEOMHT 


This routine computes the element matrices required by DERIV to evaluate the derivative 
of temperature. They include the specific heat, thermal conductivity, internal heat genera- 
tion and convection boundary condition matrices. 

It is called from MAIN. 
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GEOMHT 

















LOOK 


This is a one-dimensional interpolation routine. It is used to determine the material prop- 
erty values that are temperature dependent. These properties include specific heat, thermal 
conductivity and density and the tables are set up as properties at a specified temperature. 

It is called from MAIN and DERIV. 
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LOOK 
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MATXB 


This routine assembles the condensed global specific heat matrix from element specific heat 
matrices. 

It is called from DERIV. 
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MATXB 
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OUTPTG 


This is the routine for the optional geometry output. It is called if the user inputs KODG 
equal to 1 . 

It is called from GEOMHT. 
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OUTPTG 
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OUTPUT 


This standard routine outputs temperature, derivative of temperature, and other parameters 
at requested intervals. 

The standard routine outputs computed values in tabular form. If the user wishes, he can 
readily design an output routine conforming to the geometry of the problem (as illustrated 
for the sample problem). This is done using the variables NROWS, NCOLS and IOUT. 

OUTPUT is called from QKNINT. 


19 



OUTPUT 
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QKNINT 


This routine numerically integrates the system of ordinary differential equations written on 
nodal temperatures. It employs one of two predictor-corrector, integration algorithms. Both 
are three-stage, one-step, first order accurate explicit numerical integration procedures. 

It is called from MAIN. 
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QKNUIN 


This routine computes the new values of temperature. It performs an accuracy comparison, 
and may increase, decrease, or keep constant the step size for the next integration step. 

It is called from QKNINT. 
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qknuin 






24 



















SYMPRO 


This routine multiplies a symmetric matrix by a vector. 
It is called from DERIV. 
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III. INPUT FORMAT 


The following section illustrates the proper format and sequence for all input data. 
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IV. INPUT FOR SAMPLE PROBLEM 


This section describes the input required to generate a ten second firing of a 104 finite 
element discretization of a small axisymmetric, thrust-vector-control rocket motor as illustrated 
on the following page. The technical aspects of this test are fully discussed in Bell Aerospace 
Report 9500-920 257. Section V illustrates and provides a complete description of the output 
from COMOC using the input detailed in this section. Section VI contains the listing of the 
user-written output subroutine used for output display of this problem. 

The input is illustrated on computer forms. Where numerous cards of the same format 
are required (such as the element information cards), only two such cards are shown. The input 
is completely printed out at the beginning of output for any run by COMOC. The user should 
refer to Section III for a full description of all element and node dependent input. 
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Figure 1 . Finite Element Discretization of Thrust Vector Control Rocket Motor 
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V. OUTPUT DISPLAY AND DESCRIPTION FOR SAMPLE PROBLEM 


This section presents sample output from the rocket motor problem described in Section 
IV. A description of each portion of the output is also included. 

The first page of output is the general COMOC title page with the user supplied title 
describing the problem for which the run was made and the date. This is followed by pages 
containing a complete reproduction of input data. Included are titles, general information, element 
information, nodal information and material properties tables. This is followed by integration con- 
trol parameter information and the non-dimensionalizing constants computed in the MAIN routine. 

The following pages illustrate this information for the sample problem. 
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0.0 

0.0 

C.IOOCOE 01 

c 

5 

6 

is 

1 

C.225CCE 

C3 

0.64600E 04 

0.0 

0.0 

C.100C0E 01 

10 

19 

6 

20 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

11 

6 

7 

20 

1 

C.25CCCE 

C3 

0 • 64600 E 04 

0.0 

0.0 

0.100C0E 01 

12 

20 

7 

21 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

13 

7 

8 

21 

1 

C.35CCCE 

03 

0. 6 4600 E 04 

0.0 

0.0 

C • 1 CO COE 01 

14 

21 

8 

22 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

15 

8 

S 

22 

1 

C.4450CE 

C3 

0.64600E 04 

0.0 

0.0 

C.IOOCOE 01 

16 

22 

9 

23 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

17 

S 

10 

23 

1 

C.5CCCCE 

03 

0.646C0E 04 

0.0 

0.0 

C.IOOCOE 01 

18 

23 

10 

24 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

IS 

10 

11 

24 

1 

C.445CCE 

03 

0.64600E 04 

0.0 

0.0 

C.IOOCOE 01 

20 

24 

U 

25 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

21 

11 

12 

25 

1 

C.35CCCE 

C3 

0.64600E 04 

0.0 

0.0 

C.IOOCOE 01 

22 

25 

12 

26 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

23 

12 

13 

26 

1 

C.2500CE 

03 

0.64600E 04 

0.0 

0.0 

C.10000E 01 

24 

26 

13 

27 

1 

0.0 


0.0 

0.0 

0.0 

C.IOOCOE 01 

25 

13 

14 

27 

1 

C.2250CE 

03 

0.6 4600 E 04 

0.0 

0.0 

C.IOOCOE 01 

26 

15 

16 

2S 

1 

c.o 


0.0 

0.0 

0.0 

C.1CCC0E 01 

27 

2S 

16 

30 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

28 

16 

17 

30 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

2S 

30 

17 

31 

1 

0.0 


0.0 

0.0 

0.0 

C.IOOCOE 01 

30 

17 

18 

31 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

31 

31 

18 

32 

1 

c.o 


0.0 

0.0 

0.0 

C.IOOCOE 01 

32 

18 

IS 

32 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

32 

22 

IS 

23 

1 

c.o 


0.0 

0.0 

0.0 

C.IOOCOE 01 

34 

IS 

20 

3 3 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

35 

3 3 

20 

34 

1 

C.C 


0.0 

0.0 

0.0 

C • ICO COE 01 

36 

20 

21 

34 

1 

c.o 


0.0 

0.0 

0.0 

C.IOOCOE 01 

37 

34 

21 

35 

1 

c.o 


0.0 

0.0 

0.0 

C.IOOCOE 01 

38 

21 

22 

3 5 

1 

0.0 


0.0 

0.0 

0.0 

C.IOOCOE 01 

3 S 

35 

22 

36 

1 

c.o 


0.0 

0.0 

0.0 

C.1C0CCE 01 

40 

22 

23 

26 

1 

0.0 


0.0 

0.0 

0.0 

C.IOOCOE 01 

* 1 

36 

23 

37 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

42 

23 

24 

37 

1 

0.0 


0.0 

0.0 

0.0 

C.IOOCOE 01 

43 

37 

24 

38 

1 

c.o 


0.0 

0.0 

0.0 

C.IOOCOE 01 

44 

24 

25 

38 

1 

C.C 


0.0 

0.0 

0.0 

C.IOOCOE 01 

45 

38 

25 

39 

1 

C.C 


0.0 

0.0 

0.0 

C.10GC0E 01 

46 

25 

2 6 

3 S 

1 

c.o 


0.0 

0.0 

0.0 

C.IOOCOE 01 

4 7 

3S 

26 

40 

1 

C.C 


o 

• 

o 

o 

• 

o 

0.0 

C.IOOCOE 01 



C.IOOCOE 

01 

C * 1 OOCOE 

01 

C.1COCOE 

01 

C.IOOCCE 

01 

C . 1 OOCOE 

01 

C . 100 COE 

01 

G • 100 COE 

01 

O.IOOCOE 

01 

C « 100C0E 

01 

0 • 1 OOCOE 

01 

C . 100C0E 

01 

C.1C0C0E 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.1C0C0E 

01 

C.IOOCOE 

01 

C.1C0C0E 

01 

C.IOOCOE 

01 

0.100C0E 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCCE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

0.100C06 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

o.ioocoe 

01 

C.10000E 

01 

o.ioocoe 

01 

O.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

O.IOOCOE 

01 

O.IOOCOE 

01 

O.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCOE 

01 

C.IOOCCE 

01 

C.IOOCOE 

01 








N C T, A L DATA 


A> I S'VNNE 1PIC PF C B t E F 


FCCE 


/ ( FT. ) 


RCFT. ) 


INITIAL TEFF(F) 


I 

C.16 67CCCCF- 

-01 

C.1C822995E 00 

C.53CCCCCCE 

03 

n 

C.49S99597F- 

-01 

C.1C£22 C 55E CC 

C. 53CCCCCCE 

03 

T 

C.E2 225* 7 6 E- 

-Cl 

C.ICE22S95E CC 

C.53CCCCC0E 

03 

4 

C . 1 1 6 6 6 5 9 5 E 

CO 

C.1CF22995E 00 

C.53CCCCCCE 

03 

c 

C.1'495 9 C 9EE 

GO 

C. 916699 7 7 E- Cl 

C.53CCC0CCE 

03 

6 

C . 1 4 6 6 6 9 9 ( £ 

00 

C .E 2 3 299 76E-C.1 

C.53CCCCCCE 

03 

7 

0. 1 8 2 2 2 C C C E 

CO 

C .749999E8E-C1 

C.53CCCCCCE 

03 

£ 

0 . 1 9 9 9 9 9 9 9 F 

CO 

C .6666*9*10-01 

C . 5 3 C C C C C G.E 

03 

9 

0.216 6 69 9 EE 

GO 

C . c: 82 28999E-Cl 

C.53CCCCCCE 

03 

1C 

C. 23332595E 

CO 

C.5E229999E-C1 

C.53CCCCC0E 

0 3 

1 1 

0.25CC0CCCE 

00 

C .6666 95* 1 E — 0 1 

C.53CCCCCCE 

03 

12 

C.26666999E 

cc 

0.749/5*9 FEE - Cl 

C.53CCCCCCE 

03 

13 

C.23222996F 

cc 

C . E 2 2 2 9 9 7 6 F - C 1 

C.53CCCCCCE 

03 

14 

C.29999995E 

cc 

0.916 699 7 7E-C1 

C.53C CCCCCE 

03 

15 

C. 166 7CCCCE- 

-01 

C. 11 666995 F CC 

C. 53CCCCC0E 

03 

16 

C.*5555597E- 

-O’ 

C. 1 166 6 9 9 EE CO 

C . 5 3 C C C C C C E 

03 

17 

C. F2 229S76E- 

-01 

C . 1 1 6 6 6 9 9 5 1 C C 

C.53CCCCCCE 

03 

18 

C . 1 1 6 6 6 9 9 5 E 

cc 

C . 1 1 6 6 6 9 9 c E 0 C 

C.53CCCCC0E 

0 3 

19 

C . 1499999 EE 

CO 

C.99 C 9556**E-C1 

C. 53CCCCC CE 

03 

20 

C . 1 6 6 6 t c 9 6 E 

CO 

0. 91 669977 E- Cl 

C.53CCCCC0E 

03 

2 1 

0. 18233CCCF 

CO 

C. 8222997 6E -Cl 

C. 53CCCCCCE 

0? 

2 2 

C. 19555* 9 5E 

CO 

C. 7 45 9 <9 FEE- Cl 

C. 53CCC0C0E 

03 

2? 

0. 2 1666*9 FE 

cc 

C .6 666 9941E-C1 

C.53CCCCG0E 

03 

24 

C.23222 C 9EE 

00 

C .666699* IE- 01 

C. 53 C CCCCCE 

03 

2 5 

C. 25CCCCCCE 

0 0 

C.74999988E-C1 

C.53CCCCCCE 

03 

26 

C.266669 C 9E 

cc 

0.832299 7 6E-C.1 

C.53CCCCCCE 

03 

21 

0.28232 9 96E 

CO 

C.51669977F-C1 

C.53CCCCC0E 

03 

2P. 

0. 29999996E 

00 

C .5 C 9C5964E-01 

C. 53 C CCCCCE 

03 

29 

0.1667CCCCE- 

■01 

C . 1 3 2 2 2 9 9 9 E CC 

C. 53CCCCCCE 

03 

30 

0.499 9959 7E- 

•01 

C. 12222 9 9 9 c CC 

C.53CCCCC0E 

03 

? 1. 

0. 82 l 29916 E- 

Cl 

C.12222999E CC 

C.53CCCCC0E 

03 

3 ? 

0. U 666595E 

GO 

C.1322259 C E CO 

C. 53CCCCCC6 

03 

2 3 

C. 1499999PE 

00 

C . 1 16669. 95E CO 

C.53C CCCCCE 

03 

34 

0 . 1 6. 6 6 6 9 9 6 E 

00 

0.1CE22 C 5EE CC 

C. 5 3C CCCCCE 

03 

3 5 

0. 18332CCCE 

GO 

0 .9 9855564F-01 

C.53CCCCCCE 

03 

36 

C. 199 6 999 9 E 

00 

C.* c 98**64fc-Cl 

C. 53CCCCCGE 

03 

3? 

0 . 2 1 6 6 6 9 9 E E 

CC 

C. 95999 c 6 4 E - 0 1 

C.53CCCCCCE 

03 

3F 

C.22222995E 

CC 

C. 9999996 4fc-C 1 

C.53CCCCCCE 

03 

39 

C.25CCCCOCE 

00 

C .99959964F-C1 

C.53CCC0CCE 

03 

40 

0. 26666999E 

00 

C .89*5*9 6 4E-C1 

C. 53 C CCCCCE 

03 

41 

C. 2 E2 3 299 6 E 

GO 

C.1CF32 C 95E CC 

C. 53 C CCCCCE 

03 

42 

C. 29999995E 

00 

C.1U6 698 5E CC 

C.53CCCCCCE 

03 

4 3 

0.1667C00CE- 

01 

C. 1455959 EE CO 

C. 53 C CCCCCE 

03 

*4 

C.45595997E- 

Cl 

C.14999998E CC 

C.53CCCCCCE 

03 

45 

C .833 259 76 E- 

Cl 

C.149*9*98E cc 

0. 53 C CCCCCE 

03 

4 6 

0 . 1 1 6 6 6 9 9 5 E 

oc 

0.1 4 99999 8E CC 

C.53CCCC0CE 

03 
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47 

C. 1456656 EE 

CC 

C . 1 2 2 ? 26 6 6E 

CC 

C.53CCCCCCE 

03 

AF 

C. 166666 6 66 

00 

C.125CCCCCE 

CO 

C.53CCCCCGE 

0 3 


C.13332CCCE 

CC 

0.1 1666665E 

CO 

C.53CCCCCCE 

07 

50 

G. 19<;cq<-:cc E 

CC 

0. 116666656 

CC 

C.53CCCCCCE 

03 

5 1 

0.2 1666666 E 

CO 

C. 116666656 

CC 

C.53CCCCC0E 

03 

52 

0.22232965E 

oc 

C .1 166666 5E 

CC 

C.53CCCCCCE 

03 

5 3 

C.26CCCCCCE 

CC 

C . 1 166666 56 

CC 

C.53CCCCC0E 

03 

54 

- C.2 A 6f66 c 6E 

00 

C.l 16666656 

CC 

C.53CCCCCCF 

03 

5*5 

C. 282226 6 6E 

oc 

0 .1 25CCCCCE 

CC 

C.53CCCCCGE 

03 

5 6 

0.206966656 

CC 

C. 1222266 6E 

oc 

C. 53CCCCCCE 

03 

57 

C. 1667CCCCE 

-Cl 

C. 156? 26 666 

00 

C. 53CCCCCCE 

03 

59 

0.4966666 7E- 

-Cl 

0.156225666 

CO 

C.53CCOOCCE 

0 3 

59 

0. 8 2 2 766 7 6 E- 

-Cl 

C. 15 6 2 26666 

CC 

C. 53CCCCCCE 

03 

60 

0.116666656 

oc 

0.156326666 

OC 

C.53CCCCCGE 

03 

6 1 

C. 14666 66 EE 

CC 

C. 1416666 6 6 

CO 

C.53CCCCCCE 

03 

62 

C. 1 6 66 6 6 ^6E 

CC 

0.1 23326666 

CO 

C.53CCCCCCE 

03 

63 

G • 1 8 3 7 3CGCE 

CC 

C.125CCCCCE 

CC 

C. 53CCC0CCE 

03 

64 

0. 16966666E 

CO 

C.125CCCCCE 

CC 

C.53CCCCC0E 

03 

66 

C.21666660E 

00 

C.125CCCCCE 

CC 

C.53CCCCCCE 

03 

6 6 

G . 2 2 2 2 2 6 6 5 E 

CC 

C.125CCCCCE 

CC 

C. 53CCCCCGE 

03 

67 

C. 25CCCCCCE 

CO 

G.125CCCCCE 

CO 

C.53CCCCCCE 

03 

6? 

0. 26666666E 

CC 

C.125CCCC06 

CO 

C.53CCCCCOE 

03 

69 

0.2*3326 6 6E 

CC 

C. 12232 c 666 

CO 

C.53CCCCCCE 

03 

70 

0. 26666665E 

CC 

0.14 1 6 6 6 c 6 F 

CC 

C.53CCCCCGE 

03 
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MATERIAL FFCFERTIEE 


M/5TEP UL 

CCNC0C7 IVITY 

CCNDICT IVI TV 

SPECIFIC 

CEKSITY 

TEMPERATURE 

MJRPEF 

OIPECTION 1 

OIRECTICN 2 

HEAT 




( ETU/FT*FR*R ) 

(BTU/F T*HR*R ) 

<eio7FT#*3*R) 

( L E R/ FT ** 3 ) 

(R ) 

1 

0 . 8 3 3GC E Cl 

C.E22CCE Cl 

O . 10920 E 00 

0.50620E 03 

C* A60C0E C2 

1 

0.L667CE 02 

C.U67CF C 2 

C.13000E 00 

0.50630E 03 

0.19600E 04 

1 

0.28CCCE 02 

C.2SCCCE C 2 

C.lftCOOE 00 

0 . 506 30 E 03 

C.400C0E C A 



NCN-G I MEN 5 TONAL I 2 I NG CC N S TAM S 


7 E E R M A L CCN0LCTIV1TY = 2E.CCC 
LENGTH = 0*100 
REFERENCE TEMPERATURE* 646C.C00 
CbNSITY - 506.300 
7 ! V F FACTOR = C.C29 
SPECIFIC PEAT = C.160 


SPECIFIC FEAT MATRIX NCN-C CONSTANT = 

THERMAL CCNCLC T I V I T Y MATRIX NCN-C CONSTANT = 


INTEGRATION INFLT 


FINAL TIME(TF) = 0.C0280 

TEST CP I TER ICN (EES ) = C. 01000 

INITIAL TIME { 7 C ) = C.O 

NCMBFR CF INTEGRATION NODES* 70 

KEY-MFTFnO(KFYMTC) = 1 

MINIMUM STFP ST ZE (HMIN) = C.O 

MAXIMUM STEP SIZE(HMX) =1000.00000 

PRINT INTERVAL IDE l P ) * C.CC02E 


0.012 

0.036 
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If the user desires, he may request element output from the geometry routine (GEOMHT), 
by inputting KODG equal to one. By doing this, the user will obtain NELEM pages of geometry 
output, each of which contains the specific heat, thermal conductivity, internal heat generation 
and boundary condition matrices for a finite element. The symmetric matrices are printed in a 
lower-triangular rowwise format. 

The following sample is output for finite element number one of the rocket motor 
discretization. 
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E L E N E N T NC 


SPECIF IC 


T F E R FAl 


IMEFML 


eCLNCAFY 


ECU( *FY 


SUFFICE 


. 23 

FEAT 

C.56S1*581E-C2 C * 2 S C A3 C 6 6E- C2 C. 5<53366 10E-02 0 * 29063066 E-02 0 .2<5668307E-C2 0 . 55 33 66 10E- 02 
CCNCICT IV ITY 

0 .5C6 1 C256E Cl C.5?5e‘4E4E-C6 C.12652473E Q1 -0 . 50610266E 0 1-0 . 1 265 2483 E Cl 0.6326273CE Cl 
FEAT C E N E R A T I C N 

0.0 C.C C.O 

CCNVECT ION 

C • 26 E 4 6 524 E CC C.137E715CE CC 0 . ? 82<?<?66 5 E 00 

FLUX 

C.4C*36CF6E CC O.A2C8£6ieE CC 
CCNVECT ION 

C.O C.C C.C 0.0 0.0 0.0 


SUFF/CE FLUX 



The user may request an arbitrary number of outputs from the derivative routine (DERIV) 
by inputting LIMPT equal to the number of outputs desired (one is usually sufficient). The deriva- 
tive output details the assembly of the derivative of temperature as the program goes through the 
element loop in DERIV, and is particularly useful to ascertain that the input has been properly 
prepared. 

The DERIV output displays element node numbers and their associated temperatures, the 
thermally-averaged conductivity and specific heat, and the nondimensional element matrices and 
vectors for each element. As in GEOMHT, the symmetric element matrices are printed lower- 
triangular rowwise. Also included in this output is the step-by-step assembly of the derivative 
vector. At the completion of the element loop, the global derivative of temperature vector is 
printed out as (K)(T) + (Y) GLOBAL VECTOR.” 

Following is a sample of the derivative output for the sample problem. 


49 



FI E ► E M NUMBER 


23 


ElEf'EM FCCES 

12 13 

26 

NCCAl TEEPEE ATURES 

53C.COCC 53C.CCCC 

5 30.CCC0 

AVERAGE CCNCUCTTVITY 

= 0.871515526 Cl 


AVF P AGE SPEC I f I C E E AT 

= C.11C17C6CE CC 


SPECIFIC FEAT MATRIX 

0. 39190 1 28 E-C2 

C.2CC11845E-C2 

0.40E57159E-02 

TEERML CCNEOCT IV ITY 

M ATP IX 

C.393S9791E CO 

0. 1 5 76C02 IE 01 

C.166E6511E-C6 

Y1 VECTCR 

c.c 

C.C 

C.C 

K21 MTPIX 


0.28299665E 00 

C.26848924E CC 

C.1378715CE CC 

Y21 VECTCR . . 

0.4063fcCE6E CC 

C. 420668166 CC 

0.0 

K22 MTFIX 

C.C 

0.0 

0.0 

Y22 VECTOR 

C.C 

C.O 

c.o 

( K2 1 )*(T) VECTCR 

0. 32339149E-C l 

0.345 254C7E-CI 

c.o 

( K22 )M T ) VECTOR 

0.0 

0.0 

0*0 

(K1MITI VECTCR 

-0.11920929E-C6 

-C.im5e71E-C7 

-0.596C4645E-O7 

- ( K 1 (T ) VECTOR 

-C.33339050E-C1 

-C. 34 52535 5E -Cl 

C. 556 C46 45E-07 

-IKIIII4 (Y) VECTCR 

C.37302178E CO 

C.36633671E CC 

C. 596 04645 E-07 

-(KMT) ♦ (Y) GLOEAL 

VECTOR 

0. 1 0766CC5E 01 

0.53922063E 00 

C • 1C 1 844 C7E Cl 

0.65319 8 17E CC 

0 *101.234 1 5E 01 

C.1C225C67E 01 

0.38633371E CC 

C.O 

0.145 01 161E-07 

-0.35762787E-C6 

-0.3 57 6276 7E-C6 

C. 81 5 5638 7 E-07 

0.0 

C.556C4645E-C7 

0.0 

c.o 

C.O 

0.0 

0.0 

C.O 

C.C 

c.c 

0.0 

c.o 

0.0 

0.0 

c.o 

0.0 

c.o 

0.0 

0.0 

C.O 

0.0 

0.0 

o.c 

c.o 


0 • 200 1 1845 E-02 0 .20428 5S6E- 02 0. 40857 159E- 02 

-0 • 15760031 E 01 -0.39399821E 00 C.157C00C3E Cl 

0.0 0.0 o.c 

0.0 0.0 O.C 


0.11901960E 01 
0 .102250 77E 01 
0.1862645 IE- 06 
-0.1341 1045 E-06 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


0.99961466E 00 
0.10123434E 01 
-0.11175871 E-06 
0. 18253522 E- C6 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


0.75829774E 00 
0.E5819823E 00 
-O’. 2G489057E-06' 
-0.1152092 9E-C6 
O.C 
O.C 
0.0 
O.C 
O.C 
O.C 

o.c 



GLCEAl MATRICES 


SPECIFIC FEAT GLOE AL MATRIX - FEAT CAPACITY CCKDERSEC AT THE 


0.221311 15F- 

01 

G.66810012E- 

Cl 

0.666148556- 

01 

0.23775883E- 

01 

0.2 12765 E4E- 

Cl 

0.1855355CE- 

01 

0.25455312E- 

01 

C.1836748CE- 

Cl 

C. 53 8563 85E- 

01 

0 . 1497535 2E 

00 

C.85538447E- 

01 

C.78C05129E- 

01 

C. 10516566E 

CO 

0.66 58S654E- 

Cl 

0. 646S6829E- 

01 

0. 32045346E 

CO 

0.2128 14 16E 

CC 

C.21657170E 

00 

C. 17023331 E 

CC 

0.1625227 3E 

CC 

C. 13 774755E 

00 

C. 14562470E 

CO 

0.265 35554E 

CC 

0.26 54C464E 

00 

0. 10368866E 

CO 

0.1026 72 82E 

CC 

0. 1 02638 84E 

CO 

0. 12392223E 

CO 

0.388 14S59E- 

Cl 

0.621402 04E- 

■01 

0. 55546916E- 

01 

0.367645 CCE- 

Cl 

C.36515597E- 

01 

0.367 L3146E- 

Cl 

0. 373240636- 

01 

C.26275C72E- 

■01 

FSE GLOEAL SPECIFIC 

FEAT MATRIX 




0.45 18 524 2E 

02 

0.1456781 66 

C 2 

C.14566722E 

02 

0.4 205 23 38 E 

02 

0 .465551 30E 

C2 

0* 52645338E 

02 

0.392783516 

02 

C.54444C46E 

C2 

0. 1065CC37E 

02 

0.667583756 

01 

C.1165C649E 

C2 

C.12815012E 

02 

0.55088072E 

Cl 

0 • 1 15 48718E 

C2 

C. 11 806817E 

02 

0. 212057 78E 

Cl 

0.3 186 553 IE 

01 

C.46C88953E 

01 

0. 5874 29C5E 

01 

C.61525846E 

Cl 

0.72556550E 

01 

0. 68665662E 

01 

C.37685223E 

01 

0.376783 08E 

01 

0.5644 25 56E 

01 

C.S7356755E 

Cl 

C.57428555E 

01 

0.80655772E 

01 

0.25763229E 

02 

0. 16052636E 

02 

C. 17674C84E 

02 

0.25782453E 

C2 

G. 2 7C88806E 

02 

0.2723 81 50E 

C2 

C.26752358E 

C 2 

0. 3 6 05 88 84 E 

02 


FCDES 

0.6598 74 08E-01 
0. 18579740 E-0 1 
0. 21 52 89 89 E 00 
0.88068008 E- 01 
0.5C910864E-01 
0.13016522E CO 
0 • 12 184775 E 00 
0.2 5952506 E 00 
0. 10264564E 00 
0.93419671 E-0 I 
0 .36714818 E-01 
0.2713 37 18E-Q1 


0 .496602 80 E- Cl 
0.20451412E-01 
0.2153158eE OC 
0 • 10599726E 00 
0 • 15685028E 00 
0 • 13040 84 IE 00 
0.13017535E 00 
0.17358428E 00 
0.1026 7 3 5SE 00 
0.93442738 E-01 
0 • 36699 7 3 1 E“ 0 1 


0.26 29903 0E-C1 
0. 2254147 OE -Cl 
0 • 21 1 12550E 00 
0 • 1 14 10868E CC 
C.32C405C1E CO 
0. 15401 61 6E 00 
0.70117414E-01 
0. 10974592E CO 
C.104C9755E 00 
0. 505665236-01 
0* 36 7C72C8E-G1 


0 • 15154406 E 02 
0.53822052 E 02 
0.46448994 E 01 
0.11 354860 E 02 
0 • 19642166 E 02 
0 • 76825438 E 01 
0.82069626 E 01 
0.3653 1923E 01 
0 .97422552 E 01 
0 • 10704383 E 02 
0.27236954E 02 
0.36854507E 02 


0.20 136 810E 02 
0.48896378E 02 
0 .46443 3866 01 
0.94342051E 01 
0 • 63755 064E 01 
0.76682167E 01 
0 .768 17093E 01 
0 • 57608900E 01 
0.S7296C21E 01 
0.10701 740E 02 
0 .27 2481 54 E 02 


0 • 3802421 6E 02 
0.43589 1 72E 02 
0.47364157E Cl 
0. 8 7635746E Cl 
0. 31210489E Cl 
0.64527406E Cl 
0.14261792E 02 
C • 5 1 1 19556E 01 
C.56C63738E 01 
0.1C993054E 02 
C.27242595E 02 



The integration output is printed out at time intervals of length DELP if KPRINT = 1, 
and of slightly longer intervals if KPRINT = 0. Heading each output is the COMOC and user 
identifying titles for the problem, and the current time and integration step-size. 

The dimensional temperature and non-dimensional derivative of temperature are output 
in either tabular output form, or according to the user supplied output routine which may 
simulate the geometry of the problem. If the user does not specify a minimum step-size, the 
program will automatically compute it. Output from QKNINT will then occur when the minimum 
step-size has been determined. 

Following is the initial and final output of the sample problem of a ten-second firing of 
a small axisymmetric rocket motor. 
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C C M C C 

COMPUTATIONAL CONTINUUM MECHAN ECS 
TFEPMAL ANALYSIS VARIANT 

STANCARD TEST CASE FCR CCMCC 

104 FINITE ELEMENT C I SCR ET I Z A T ICN CF A SMALL THRLST VECTCR CCNTFCL 
RCCKET MOTOR, LSING TRIANGULAR AXISYMMETPIC RING ELEMENTS WITH VARIABLE 
CCNVECTIVE 8CUNCARY CGNCITIONS APPLIEC TC INTERIOR SURFACE • TPAMSIENT 
BEFAVICR CURING A 1C SECOND FIRING DETERMINE!) LSING TEMPERATURE 
CFFENDENT MATERIAL PROPERTIES AND CKNINT ALTCMA TIC INTEGRATOR* 


tn 

co 


TIME = 0.0 

HOURS 

CURRENT INTEGRATION STEP SIZE = 

0.10000 E -05 FOURS 



X1-CCCRCINATE( FT.) 
0.02 0.C5 

0.C8 

C.12 C.15 C.18 0.20 

0.22 0.23 0.25 

0.27 

C. 30 


TEMPERATURE! DEGREES RANK1NE) 











530.00 

530. CO 

530. CC 

530. CO 











53 C .GO 

530.00 

530. CC 

53C.CC 

53C.CC 









53C.00 

530.00 

530. CO 

530. CO 

530.00 

53C.CC 

5 20.CC 







530. CO 

53C.CC 






530. CO 

53 C. 00 

530.00 

530.00 

530.00 

530.00 

530.00 

530.00 


530. CO 

530.00 

530.00 

530 . OC 

530. CC 


530.00 

530.00 

530.00 

530.00 

530.00 

530.00 


53C.C0 

530. CO 

530. CO 

530. CG 

530.00 

53C.CC 

530. CC 

53C.C0 

53G.00 

530.00 

530.00 

530.00 

530.00 

530.00 

53C.C0 





53C.CC 

530. CC 







530. CC 

530. CC 






530. CO 

530.00 





530.00 

530. CO 








53C.C0 

530.00 



530.00 

530.00 










530.00 

530.00 

530. CO 

530.00 





530.00 530.00 


TIME CEPJVATIVE 

OF TEMP 

ERA TORE ( NGN 

-DIMENSIONAL) 


0.00 

0.00 

-O.CC 

-0.00 




-o.co 

0.00 

O.CC 

c.oc 

-c.cc 



0.00 

0.00 

0.00 

-C.00 

c.cc 

C.CC 

C.CC 

-0.00 

0.00 

C.00 

O.CC 

-O.CO 

-O.CO 


0.00 

0.54 

1.08 

i.ce 

1.15 

-C.CC 

o.co 

C. GO 





l.CC 

-O.CC 







0.76 

-C.00 







C.86 


CKNINT 

CUMULATIVE NUMBER OF TIMES STEP SI2E WAS DECREASED = 


-0.00 

-0.00 

0.00 

-0.00 

-0.00 

0.00 

-0.00 

-C.CC 

0.00 

0.00 

0.0 

0.00 

0.00 

0.00 


0.00 

0.00 

0.0 

0.00 

0.00 

0.00 

-O.CC 

c.cc 

-0.00 

1.01 

0 

0.00 

1.02 

AND 

-O.QO 

1.02 

INCREASED 

0.00 

1.01 

0.00 

0.86 

0 

O.CC 

0.76 

C. 38 


QKMNT 

TIME AT WHICH ACCLRACY TEST WAS FIRST FA I LED = 
STEP-SIZE AT PREVIOUS TIME STEP 
NEW MINIMUM STEP-SIZE ( FM IN ) 


C.845S59E-04 HOURS 
Q.16CC00E-04 HOURS 
C.16CCC0E-04 HCUPS 



COMPUTATIONAL CCNTINLUM MECHANICS 
THERMAL ANALYSIS VARIANT 


STANCARO TEST CASE FCR CCMCC 

10A FINITE ELEMENT C I SCR ET IZA T ICN CF A SMALL THRLST VECTOR CCNTFCL 
PCCKET MCTCR, USING TRIANGULAR AXISYMMETRIC PING ELEMENTS VITH VARIABLE 
CONVECTIVE BCUNCARY CCNCITIONS APPLIED TO INTERIOR SURFACE. TRANSIENT 
BEHAVI CF DURING A 1C SECOND FIRING DETERMINED L SING TEMPERATURE 
CEFENDENT MATERIAL PROPERTIES AND CKNINT AUTOMATIC INTEGRATCR. 


TIME = 0 .0028 E S FOURS CLRRENT INTEGRATION STEP 


X1-CCCRCINA1E(FT 

.1 





0.02 

0.05 

0.C8 

0.12 

C.15 


0.18 

TEMPERATURE! CEC-REES RANKINE) 




686.92 

656.52 

651.54 

672.16 




727.01 

728.82 

723.38 

7C4.75 

791. C7 



1009.22 

587.57 

577.19 

546.25 

£24. CC 

841.38 

852.22 

657.43 

1541.16 

1505.54 

1454.20 

1456.06 

1124.46 


555.34 

1873.09 

1837.69 

1828.46 

1757.59 

me. 2i 

1241.58 

1367.61 





2C87.62 

1502. CC 
2218. e9 

2123.66 

2565.54 


TIKE CEF IVA1IVE OF 

TEMPERATURE ( NON-O IMEN SICNAU 


0.04 

0.07 

0.06 

0 .06 




0.11 

0.20 

0.19 

0.18 

c.c< 



in 

H 

• 

O 

0.21 

0.30 

C.28 

o.ie 

0.05 

0.14 

0. 05 

0.09 

0.21 

0.21 

0.2 0 

0.26 


0.14 

0.02 

0.06 

0.C6 

C.C6 

o.ie 

0.18 

C. 20 





0.C6 

0.12 

0.03 

C. 12 
C. 04 

CKMNT 
CUKUl^T 1VE 

NUMBER 

OF TIMES 

STEP 

SIZE WAS 

CECREASEC 



SIZE = 0.30A00E-O3 FOURS 


0.20 

0.22 

0.23 

0.25 

0.27 


0.30 







694.35 






778.56 

735.89 

810.51 

756.69 

755.38 

799.58 

854.22 

818.04 


858.12 

799.90 

799.10 

848.60 

919.50 

1168.99 

970.45 

1229.00 

1145.33 

1147.29 

1227.54 

1345.12 

1854.71 

1561.37 

1996.51 





2105.65 

2289.90 


2417.28 



2385.04 

2571.45 



2925.73 

2730.39 

2712.61 

2906.36 





3264.56 

3257.17 



0.03 

0.02 

0.03 

0.05 

0.04 

0.04 

0.04 

0.04 

0.13 


0.13 

0.13 

0.13 

0.13 

0.13 

0.17 

0* 08 

0.26 

0.29 

0.27 

0.22 

0.18 

0.11 

C.C6 

0.02 


4 ' 



0.14 

0.03 


0.16 



0.20 

0.03 



0.04 

0.22 

0.24 

0.04 





0.04 

0.04 





1 

ARC 

INCREASED = 

22 





VI. USER-WRITTEN OUTPUT SUBROUTINE FOR SAMPLE PROBLEM 


Contained in this section is a FORTRAN listing of the user generated subroutine which 
presents output in a format geometrically similar to the sample problem, as was illustrated in 
Section V. To operate this output feature for any problem, the user sets NOUTP to any non-zero 
integer and provides COMOC with a routine titled SUBROUTINE OUTPUT. The variables 
NCOLS and NROWS are then used to set up a specific output format. 
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CJ1 

CD 


FCPTPfN TV G LEVEL 19 


CITPU C*TE = 72300 13/ 12/46 


00 01 

c 

SUBROL TINE GLTPLT 


c 

Q 

LSER GENERATED OLTPLT RCLTINE FCR RCCKET MCTCF 

CCC2 


COMMON / CKNI / VV (2 *7CG) * ZZ(2,7C0), XX(25) 

GO C3 


C0MMCN/CCMR/CS(1CC0,6),AKS(1CCC,6) ,FGl (1000 ,3) ,ZM< 1000,3 ) , 

1 AK2S(1CCC,6) ,\M2(1CCC,3) ,YM( 1000 ,2 ) ,FC(1000), C< 1000 ) , 

2 TR( 1CCC) * IK ( 1 C CC ) * F ( 1 C 00 ) »RCCCP(700) ,ZCCCR ( 700 ) * T( 700), 

3 C(7CC)»TTABC(2C*6) ,T1ABK(20*6),TTAERC(20,6), TFIX(IOO), 

A CTAB(2C,6) »AKTAB(2C,6) , PCI A E (20 ,6 ) ,F EAC ( 6 , 18 ) , 

5 AlC ,AKC ,TRC ,TCCNST,F1IN .AK2TA8 (20 ,fc) 

0004 


COMMON/ CO MI/MELEM(1CCC) ,KBND I 1CC0I , I NOCE (3000 ) ,NT A eC <6 ) , NT ABM 6 ) , 

1 NTABRC(t) *1 CL T (2C,30) ,NCCLS(30) , NR CWS , I F I X ( 100 ) , NF IX, 

2 NTITLE,NELEM,NNCCE,NCCCPC,NN,NNS,IFIRCI ,KCCA , KCCC , L IMPT 

CCC5 


CIMENSICN Z ( 7CC ) 

CCC6 


CATA KCLNT/O/ 

0007 


EQUIVALENCE (XX(24),N),(XX(21),NCE) 

ocoe 

1 

FORMAT ( 1 HC / / 1 TIME =*,F12.6,' HCURS* , 

1 BX, 'CURRENT INTEGRATION STEP SIZE =*,E15.5,' HCURS') 

occs 

1C 58 

FORMAT! 1H, A0X,'C 0 M C C ' , //30X , 'CCMPL T AT I C N A l CONTINUUM MECHANIC 
IS ', /3AX, 'THERMAL ANALVSIS VARIANT',//) 

0010 

1055 

FORMAT ( IE 1 1 8A4 ) 

0011 

1100 

FORMAT (IF ,15(F7*2,1X)) 

0012 

1101 

FGRMAT(1F0,//,1H , 'TE MPERA TURE (DEGREE S FANKINE)') 

CO 13 

1103 

FCRMATt 1H0,'TIME DERIVATIVE CF TEMPER A TORE t NC N-C I M ENS ICN AL ) • ) 

0014 

1120 

FORMAT ( 1H0, • Xl-OCORD I NATE(FT.) • ,/IH , 

1 5(F?.2,1X),BX,6{F7.2,1)),8X,F7.2) 

0015 

1200 

FORMAT ( IP ,32X»F7.2*65X*F7.2) 

0016 

1 2C5 

FORMAT ( IF ,6(F7.2,1X),A8X,2(F7.2 ,1X)) 

0017 

1210 

FORMAT ( IF , FOX, 11 (F 7. 2, IX)) 

0018 

1220 

FC B MAT ( IF »5(F7.2»1X) ,EX,6(F7.2,1>) ,8 X , F 7.2 ) 

0019 

1230 

FORMAT ( IF ,A(F7.2,1X) ,EX,F7.2,49X,F7.2) 

0020 

1240 

FORMA T( IF ,32X,F7.2,«X,6<F7.2,1X) ,8X,F7.2) 

0021 

1250 

FORMAT ( IF ,32X,2IF7.2,1X) ,A8X ,2 < F7. 2 , 1 X ) ) 

0022 

1260 

FORMAT ( IF ,F0X,2(F7.2,1X) ,32 X ,2 (F7. 2 , 1 X ) ) 

0023 

1270 

FORMAT ( IF ,A6X,2(F7.2 ,1X) ,16X ,2 ( F7.2 , 1 X ) ) 

0024 

1280 

FCRMA T( IF ,56X,4(F?.2,1X)) 

0025 

1250 

FORMAT ( IF ,6AX,2IF7.2,1X)) 

0026 


WR)TE(6*105€) 

0027 


CO 10 t = l ,NT ITLE 

0028 

10 

WR ITE (6,1059) (HEAD ( I ,J) ,J=1,IE) 

0029 


TME = XX( 15) * TCCNST 

0030 


SS = XX(1A) * TCCNST 

0021 


WRITE (6,1) TME , S S 

0032 


CO 5 1*1, NOE 

0033 

c 

Z ( I) =YY(N,t)*C(I)*TRC 

0034 


00 2 1=1 ,NNODE 

0025 

2 

Z COCR ( I ) = ZCCCR ( I ) * ALC 

0036 


N 1 = NCCL S ( 6 ) 

0037 


WRITE(6,1120)(ZCCCR(ICLT(6 ,J)) ,J=1,N1) 

0C3 8 


WRnE(6,UCl) 

0029 


CO 301 I=1,NPCV$ 

0040 


N 1 = NCCL S ( I ) 

0041 


GO TO ( 11,11,13,1*, 15,16*1 7, 18, IS, 20, 21 ,22,23) ,1 

0042 

11 

WR ITEI6, 1100)<Z( IOL T <1 ,J) ) ,J = 1 ,N1 ) 

0043 


GO TO 301 

0044 

13 

WR !TE(6,1200)(Z(ICLT(I ,J) ) 

0045 


GO TO 301 


90003 


10009 


PAGE 


0001 



FCPTP^ IV G LEVEL 


CLTPL 1 


C*TE 


72300 


0C46 
CC47 
C C 48 
CC49 
CC 50 
0 C 5 1 
00 5 2 
C Q 5 3 
CC 54 
CC 55 
C C 56 
CC 57 
CC58 
C C 59 
CCfcO 
C C 6 1 
CC62 
CC6 3 
CC64 
0C6 5 
C 066 
CC67 
CC 68 
C C 69 
C C 70 
C C 7 1 
CC72 
CC 73 
CC74 
CC75 
CC 76 
CC77 
C C 78 
CG79 

cceo 

C C 8 1 
0C62 
CC 63 
C C 64 
0 C 6 5 
CC 86 
C C 67 
CC 68 
C C 89 
CC90 
0C91 
C C 52 
C05 3 
C054 
CC55 
CC 5 6 
C C 5 7 
C C 58 
C059 


15 


14 M*ITE<6 V 1205HZ(IQLTU ,*)> ,J = 1 ,M) 

GO TO ?C1 

15 WP I T E ( 6 , 1210)(Z(ICLT(J , J } ) , J = 1 ,M ) 

GO TO 301 

16 WR ITE(6, I22G) ( Z( IOLT (I , J) ) ,J = 1 ,M) 

GO TO 30 1 

17 WP ITEI6 ,1220 H Z( IOLT (I , J) ) ,J = 1 ,M ) 

GO TO 301 

18 WP l TE ( 6 , 1240 > < Z( IOLT (1 , J ) ) ,J = 1 ,M) 

GO TO 301 

19 WR ITE< 6 , 1250MZUOLTU ,J)) , J = 1 ,M ) 

GO TO 301 

20 WR ITE( 6, 126C HZ (IOLT (I ,J)) ,J = 1 ,N1) 

GO 70 301 

21 ViPITE(6»127C)miCLT(I f J)) V J«1 ,m 
GO TO 3 C 1 

22 WRlTE(6fl2£CHZ( IOLT (I , J ) ) , J*1 ,M) 

GO TO 201 

23 WR ITE(6, 1250 H Z( IOlT(I , J) ) , J = 1 ,M ) 

301 CONTINUE 

WR I T E ( 6 » 1 103 ) 

CO 211 I=1,NRCW$ 

N 1 = NC CL S ( I ) 

GO TC ( 21 , 31 , 32, 24,35, 36, 37, 36 , 39 , 40 , 41 , 42 , 43), I 
31 WRITE* 6, 1100) ( ZZ ( K , ICt 7 (I , J > ) , J«1 , M ) 

GO TC 211 

23 WR !TF(6, 12C0)(ZZ(N,ICLT(I , J) ) , J = l , M ) 

GO TO 311 

34 WRITE(6,12C5MZZ(N,ICLT(I,J) ) » J= 1 , M ) 

GO 70 311 

35 WR !TE( 6, 12 10 HZ Z(K,irLT (I ,J) ) ,J = 1 ,M) 

GO TO 211 

36 WR lTE(fc , 12 20) 122 IN ,ICLT (I , J) ) , J = 1 ,M ) 

GO TO 311 

3 7 WP ITE(6,1230)(Z2(N,ICLT (I , „) ) ,J = 1 , M ) 

GO TO 211 

3 8 WR !TE<6, 1240HZZ(N,I0LT(I , J) ) , J = 1 ,N1 ) 

GO TO 211 

39 WR ITE(6,1250MZZ(N,I017(I , J ) ) , J=1 ,N1 ) 

GO TO 2 11 

40 WR IT EC 6, 126 0 HZ 2(N,ICL7<I ,J) ) ,J=1 ,M) 

GO TO 211 

41 WP ITE(6, 127C )(ZZ (NtlCLT (I ,J) ) , J=1 ,N1 ) 

GO TO 211 

4 2 WR ITE ( 6 * 12 8C H Z Z ( N , I Ot l.( I , J) ) , J=1 , M ) 

GO TO 211 

4 3 WR ITfc ( 6, 1 25C HZZ (N , IOLT (I, .j ) ) ,J*1,M) 

211 CONTINUE 

CO 2 2 C 1 = 1, NNCCE 
32C ZCnCK(l) = ZCCCR(I) / ALC 
K0UM = KCLM+1 

IF <LI*PT .CT* KCLNT) KCC4 = 1 

RETLPN 

ENC 



